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I. INTRODUCTION 

In introductory courses on cosmology, it is common 
to study the evolution of the scale factor assuming 
that the universe is filled with a perfect fluid with 
equation of state 

p^wp, (1.1) 

where p is the pressure, p is the energy density (in- 
cluding, of course, the density of mass multiplied by 
c^) and w is a dimensionless constant. This covers, in 
particular, the very important cases of dust {w — 0) 
and radiation (w = 1/3). 

Although it is relatively easy to write down the rel- 
evant equations (see Sec. HIl below), the time evolution 
is usually discussed only for a few specific cases. Of- 
ten, textbooks present the exact solutions correspond- 
ing to dust and radiation, and sometimes even stiff 
matter (w = 1), in the absence of a cosmological con- 
stant A. This, however, leads some students to won- 
der how general are the results so obtained. In the 
literature, one can also find analytic solutions for an 
arbitrary value of w (see, e.g., Refs. [H-d), but a de- 
tailed analysis of the properties of such solutions is 
somewhat lengthy and time-consuming — certainly 
not appropriate for introductory lectures. Further- 
more, it must be noted that it is not terribly inter- 
esting to know the analytic solutions for simple cases 
with a definite w, given that none of them applies 
to the entire history of the universe, and that an ex- 
act solution in which the universe in different eras is 
dominated by different components is definitely out of 
reach. 



Even leaving such considerations aside, the exact 
solutions that one can find in the literature usually 
do not cover the possibilityof a non- vanishing cosmo- 
logical constant (see Refs. '3-'6' for notable exceptions, 
which are, however, restricted to the case of dust). 
Yet the current observational data not only do not 
exclude, but appear instead to support a cosmologi- 
cal constant''"- (although other possibilities, involving 
the broad concept of "dark energy", are often con- 
siderediSiii) . Apart from very special cases, finding 
exact solutions even for a specific value of w and a 
non-zero A is usually a mathematical problem well 
beyond the capability of most of the students, and 
would anyway take time that is better used to present 
other topics of greater physical interest. 

One of us was facing these problems while teach- 
ing a course whose aim was to give a broad, yet ac- 
curate presentation, in 28 hours, of modern physi- 
cal cosmology to undergraduate students majoring in 
mathematics, computer science, and engineering. His 
strategy was to entirely avoid exact solutions, replac- 
ing them by a qualitative study. Remarkably, such 
qualitative analysis, together with the asymptotic be- 
haviour of the scale factor in regimes where it is very 
easy to solve the equations, allows one to gather most 
of the relevant information without dealing with un- 
necessary mathematical technicalities. Although it is 
traditional, in elementary courses, to cover only the 
range ^ w ^ 1, this method can be easily applied 
for all values of w. A paper published a few years 
agoi^ contains a nice discussion of the physical prop- 
erties of matter, and of their cosmological relevance 
and implications, with w in different ranges. Here, we 
present a complementary analysis, showing how such 
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kinds of matter affect the evolution of tlie universe at 
large scales. 

Qualitative techniques are well-known and em- 
ployed in cosmology. ^"^^^^ However, no elementary dis- 
cussion, suitable for being used in an introductory 
course, is available to our knowledgcj^ In this pa- 
per we give a very simple presentation, based on the 
often noted analogy between the scale factor and the 
coordinate of a one-dimensional mechanical system. 
Surprisingly, we could not find such a treatment in 
a fairly large sample of books on cosmology that we 
have analysed, so we believe it is worth bringing it to 
the attention of the community of physics teachers. 

In Sec. In] we write down the basic equations. The 
mechanical analogue, and an outline of the qualitative 
analysis, are discussed in Sec. IIIIl while Sec. |IV] con- 
tains the detailed study for physically interesting com- 
binations of w and A. In Sec.|V]we present some final 
comments and possible extensions of our approach. 



II. EQUATIONS 

Our notations and settings are standard. We 
consider a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) model, whose t = const spatial sections are 
homogeneous and isotropic. The spacetime metric is 



dr^ 



sm 



1 — kr"^ 

(2.1) 

where a{t) > is the scale factor, k = —1, 0, or 
1 according to whether the spatial geometry is hy- 
perbolic, flat, or spherical, and r, 9, ip are comoving 
coordinates. Note that r is dimensionless. We shall 
restrict our analysis to the case in which the energy 
density p is strictly positive; an extension to p ^ is 
mathematically straightforward. 

The equations describing the evolution of this model 
in the presence of a perfect fluid and a cosmological 
constant A are: 

• The Friedmann equation 



• A barotropic ° equation of state 



(2.4) 



usually chosen of the type (II. ip . 

Overall, we have three equations for the three un- 
known functions a{t), p{t), and p(t). Replacing the 
equation of state into Eq. (|2.3p one finds 



p a 
which can be immediately integrated to obtain 



P ■ 



3(l+iu) 

Po a-0 

„3(l+u.) 



(2.5) 



(2.6) 



where the subscript corresponds to some arbitrary 
reference time to- This can be replaced into the Fried- 
mann equation (j2.2l) to get a single differential equa- 
tion for the unknown function a{t): 



■2 /I /-< 3{l+w) , 

a AnGpo Oq 1 



3c2 



ai-.-— (2-7) 



where we have multiplied by 0^/2 and rearranged the 
terms for later convenience. 



III. EQUIVALENT MECHANICAL SYSTEM 

Equation ()2.7|) is formally equivalen t "'^^1^° to 
a2 



+ Via) = E, 



(3.1) 



expressing the conservation of energy for a particle 
with unit mass in one dimension, with coordinate a 
and total energy 

E = ~kc^/2, (3.2) 
subject to forces with potential energy 

A 



Via) = -- 



+3to 



(3.3) 



SttG 



where, as usual, a dot denotes the derivative with re- 
spect to t\ 

• The equation for the conservation of energy during 
the cosmological expansion, 



p-|-3-(p + p) = 0; 
a 



(2.3) 



where 



3c^ 



B := - — . (3.4) 



The two contributions to V are plotted in Figs. [T] 
and[21 One notes that qualitative differences arise for 
the critical values w = —1/3, w = —2/3, and A = 0. 
Another exceptional value is w = —1, because in this 
case also the first contribution becomes quadratic, so 
V{a) = {B- A)a'^, with A = AirGpo/ {ic^), and the 
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FIG. 1: The contribution ~A/a^+^'^ to the potential en- 
ergy for different values of w. 
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FIG. 2: The contribution Ba to the potential energy for 
different values of A. 



evolution is entirely controlled by an effective cosmo- 
logical constant 



Ae = A + 6A/c2 = A + SvrGpo/c^ 



(3.5) 



Given V{a) and a value of E, Eq. (|3.ip defines a 
set of points in the phase plane (a, d). This set corre- 
sponds to one or more phase curves, each associated 
with a possible solution of the equation of motion with 
energy E^^^ Since the kinetic term d^/2 in Eq. (13. ip 
is very simple, such curves can be easily drawn know- 
ing the plot of V{a). In particular, it is worth noticing 
that the Friedmann equation (|2.7p . hence Eq. (|3.ip . 
is invariant under time-reversal, so all the diagrams 
in phase space will be symmetric with respect to the 
horizontal axis. As is well-known from classical me- 
chanics, from the knowledge of the phase curves one 
can get some useful information about the dynamical 
behaviour of the system. 



There is only one place where the qualitative study 
applied to the cosmological case differs from the anal- 
ogous discussion in particle dynamics. Usually, for 
a point particle, the initial conditions affect only the 
value of the total energy E but not the potential en- 
ergy, which is fixed. Initial conditions associated with 
different values of E will then produce different phase 
curves. For a given V{a)^ there will be infinitely many 
of them. On the contrary, in the cosmological case 
there are only three possible values for E (c^/2, 0, or 
— according to the spatial geometry), so if V{a) 
is fixed there is only a finite number of phase curves. 
However, A depends on po and oq, so it can also vary 
from solution to solution, which makes V{a) depen- 
dent on the initial conditions. Therefore, in order to 
draw a representative sample of phase curves, one may 
need to consider, for each value of E among the three 
possibilities, also different values of A which can lead 
to different qualitative behaviours. 

From the potential energy V{a) one can also in- 
fer the behaviour of another quantity of cosmological 
interest: The second time derivative d of the scale 
factor, which is related to the deceleration parameter 
q ~ ^a'djc?. Indeed, taking the time derivative of 
Eq. dSj]) one finds 



dF 
da ' 



(3.6) 



which is, of course, Newton's law of motion. With 
V(a) given by Eqs. p.3p - p.4p . this is equivalent to 
Friedmaim's equation in the form 



47rG 



Ac^ 

"iw) a H — — a. 



(3.7) 



where p is given by Eq. (12. 6p . By the assumption 
p > 0, the first contribution on the right-hand side of 
Eq. p.7p is strictly negative if and only if w > —1/3, 
and vanishes for w = —1/3. Hence, a positive d, as 
currently observedj^i^i^ requires a positive cosmolog- 
ical constant or exotic matter with w < —1/3. This 
supports a study of FLRW models with a wide range 
of combinations of w and A. 



IV. QUALITATIVE STUDY 

Given the existence of the critical values w = —1/3, 
w — —2/3, w = —1, and A = 0, there are 21 possi- 
ble classes of models. We shall group them together 
according to the value of the parameter w. How- 
ever, only the physically most interesting cases will 
be treated in detail. We do not discuss at length 
the physical meaning of the different possibilities; see 
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Ref. [T^ for an instructive presentation of the physics 
underlying different models and for an extensive list 
of references. 

Let us first explain the conventions used in the fig- 
ures. For each class of models considered, there are 
three interesting plots to show: One for the effective 
potential energy V{a), one for the phase portrait (the 
set of phase curves for a given V), and one for the 
corresponding time evolutions of the scale factor a{t). 

The potential energy is always drawn above the cor- 
responding phase portrait, and the units along the a- 
axis are the same for both diagrams. We omit the ar- 
rows that are sometimes drawn on the phase curves, 
which point to the right when d > (upper half of 
the phase plane) and to the left when d < (lower 
half). The values k = —1 and k = 1 correspond to 
the values and ~(? jl of the total energy, and to 
the upper and lower horizontal lines in the diagram 
for V{a\ whereas fc = corresponds to = and 
to the a-axis. For convenience, we associate different 
colours with different values of k: Red, black, and blue 
for A; = — 1, fc = 0, and fc = 1, respectively. All the 
curves pertaining to a value of fc are always drawn in 
the corresponding colour. We use colours for greater 
clarity, but all the information can be retrieved also in 
grayscale. A comparison between the phase diagram 
and the corresponding curve for V{a) makes it very 
easy to understand which phase curves correspond to 
a given value of the energy E\ therefore, we do not 
introduce any label for this purpose. Sometimes, it is 
useful to draw plots for V(a) with different values of 
the parameter A (see, e.g.. Fig. 15]). In this case we use 
different styles (solid, dotted, dashed) both for such 
plots and for the corresponding phase curves. 

The colour and style of the curves in the diagrams 
for a(i) are the same as for the corresponding phase 
curves. Since here the correspondence with the value 
of E is less obvious in greyscale, we have added a few 
labels. Whenever all the three curves ait) correspond- 
ing to the three possibilities fc = — 1, fc = 0, and fc = 1 
are drawn for the same t^(a), they are identified by a 
"— " , no label, and a , respectively. From the di- 
agrams below one can see that this labelling involves 
only solid curves. Whenever there are two discon- 
nected curves that are related by time-reversal, corre- 
sponding to two alternative histories of the universe 
with opposite signs for d (as it happens, for exam- 
ple, for w > —1/3 and A = in the cases fc = — 1 
and fc = 0, which admit solutions that contract to a 
big crunch together with those that expand from a 
big bang, see Sec. IIV A II below), only the one with 
d > will be plotted. The other (physically irrele- 
vant) curve a{t) with d < can be obtained simply 
by reflection through a vertical axis. Also, note that 



the time t = 0, as well as the value of t correspond- 
ing to the big bang, has nothing special, because the 
Friedmann equation (12. 7p and Eq. (jS.ip are invariant 
under time translations. This allows us to plot the 
curves starting from the same point (as in Fig. H]) , but 
when this is done, it is only for graphical convenience. 
For the same reason, possible intersections between 
curves have no physical meaning. 

A. w > -1/3 

This applies to all cases of "ordinary" matter, in- 
cluding in particular dust (w = 0) and radiation 
{w = 1/3), but also to several exotic substances, like 
e.g. stiff matter {w = 1). Since l-|-3it; > 0, the poten- 
tial energy —A/a^~^^^ is singular at a = 0, increases 
monotonically with a, and tends to for a — ^ -l-oo. 
In all these models there is a big bang and/or a big 
crunch, i.e., a = at some time t~ in the past and/or 
i+ in the future. For a — >■ 0, the terms in A and in fc 
are negligible in Eq. (12. 7p . which becomes asymptoti- 
cally 

• 2 A 3(1+11)) , 

2 ^ 3c2 ai+3«' ■ ^ ■ ' 

This is immediately integrated to obtain 

(4,2) 

valid for t close to t^. The negative signs corresponds 
to the big bang, the positive signs to the big crunch. 
As the big bang/big crunch are approached, d — > ±oo. 

1. A = 

The potential energy and the phase portrait are 
shown in Fig. |31 For fc = 1 the scale factor a{t) starts 
from zero with a big bang, reaches a maximum, and 
then symmetrically recollapses to zero in a big crunch. 
For fc = Oorfc = — 1, the expansion continues forever 
but slows down, and for a — > +oo, d ^ — fcc. There 
is also a physically irrelevant time-reversed process in 
which d < 0. Interestingly, the diagram for the three 
possible behaviours of a(t), see Fig. SI is qualitatively 
indistinguishable from those presented in many text- 
books for the case w = 0, although here the value of 
w > — 1/3 is arbitrary. 

The curves corresponding to fc = in the phase 
portrait of Fig. [3] and in Fig. 2] mark the boundary 
between two qualitatively different regimes — one in 
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FIG. 3: Potential energy and phase portrait for u; > —1/3, 
A = 0. 



a 




FIG. 4: Time dependence of the scale factor corresponding 
to the phase curves of Fig. O 



2. A < 

A negative cosmological constant corresponds, in 
Newtonian language, to an attractive universal force. 
The potential energy is the one plotted in the upper 
part of Fig. [51 and leads to the phase portrait shown 
in the lower part. The phase of expansion is always 
followed by a recoUapse and a big crunch (so the plots 
of a{t) all look, qualitatively, like the one in Fig. |4] 
that is labeled by "+"), independently of the value of 
k. The asymptotic behaviour near the big bang and 
the big crunch is still given by Eq. (|4.2|) . 



3. A > 



which models recollapse and one in which models ex- 
pand forever. In general, qualitatively different be- 
haviours occur when the potential energy has maxima 
(even improperly, when it is attained for a — > -l~oo, as 
in Fig. [3]). We shall call separatrix the set of points 
in the phase space corresponding to a total energy 
equal to such a maximum value of V . This set al- 
ways separates regions of the phase space containing 
phase curves corresponding to qualitatively different 
behaviours. 



The potential energy is plotted in the upper part of 
Fig. [S] It is always negative, and for a — a, with 



(1 -H3u;)A\ 3(1+™) 



ao 



-2B 



Ac4 

attains its maximum value 



(4.3) 



3(1 + w) 2 / -2B 
V{a) = ^-^ ^3(1+™) 



2 Vl + 3w, 

a^Ac^ (4.4) 



2(1 + 3w) 
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FIG. 5: Potential energy and phase portrait for w > —1/3, 
A < 0. In all cases, the time dependence of the scale 
parameter is similar to the one labeled by a "+" in Fig. |4] 



The phase portrait is shown in the lower part of Fig.|6l 
For fc = — 1 and fc = the universe starts with a big 
bang, with asymptotic behaviour given by ()4.2|) . and 
traverses a phase in which a decreases, followed by one 
in which it increases. Asymptotically, for large values 
of a, Eq. p.7|) becomes 



Ac^ , 
2 6 

and gives the de Sitter expansion 



(4.5) 



(4.6) 



where if is a constant (see the two leftmost solid 
curves in Fig. [7]). Of course, there are also the time- 



FIG. 6: Potential energy and phase portrait for w > —1/3, 
A > , and three possible values of A. The phase curves 
corresponding to fc = — 1 and fc = are drawn only for 
the potential energy represented by the solid curve, since 
in the other two cases they look qualitatively the same. 



reversed solutions (not shown in Fig. [7]) in which the 
universe contracts from an asymptotic de Sitter phase 
at t — > — oo and ends with a big crunch. 

For fc = 1 the situation is more complex, because 
the behaviour of the scale factor depends on the spe- 
cific values of the parameters. There are three qual- 
itatively different possibilities, according to whether 
the maximum value V{a) of the potential energy is 
greater, equal, or smaller than — c^/2. 

If V{a) > —(? jl (which is, of course, an algebraic 
condition involving ap, po, and A, and corresponds to 
the dashed curves in Fig. 15]), there are two types of 
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FIG. 7: Time dependence of the scale factor corresponding 
to the phase curves of Fig. [G] 



behaviour, according to the initial conditions. Either 
the universe starts with a big bang, reaches a maxi- 
mum size (smaller than a), and then recollapses in a 
big crunch, as in the models with A ^ 0, or it starts at 
t — ^ — oo with an asymptotic de Sitter phase in which 
a decreases as 



(4.7) 



reaches a state where a has a minimum value (greater 
than a), and then bounces towards another asymp- 
totic de Sitter phase for t -l-oo, in which a{t) in- 
creases according to (|4.6I) . In Fig. [7] these time depen- 
dences are represented by the two dashed curves. 

For V{a) = — (dotted curves in Fig. ^ there 
is a separatrix in the phase plane. This is the union 
of five phase curves, one of which is degenerate (i.e., 
reduces to a single point; see the dot in Fig. [6]) and 
corresponds to an unstable static model with a{t) = 
a = ao- From Eq. (|4.3p . this happens when 



The other phase curves pertaining to this case corre- 
spond to models that start with a big bang, or with an 
asymptotic collapsing de Sitter phase, and end up ap- 
proaching asymptotically the previous static universe 
for t — > +0O, or with their time-reversed version. In 
particular, the case in which a model that is static in 
the asymptotic past evolves towards an expanding de 
Sitter universe was discovered by Lemaitrer^ and sup- 
ported by Eddington.^^'^^ (See also Ref.^for a mod- 
ern assessment of the stability of the static model.) 
In Fig. [7] these time dependences of the scale factor 
are represented by the dotted curves. As already an- 
nounced, only the curves with d ^ are shown. 

Finally, if V^(a) < — c^/2 (solid curves in Fig. |6l 
and curve labelled by a in Fig. [7]), we have a 
situation qualitatively similar to the one discussed in 
the cases fc = — 1 and fc = 0. The only noteworthy 
difference is that now it is possible for V{a) to be 
arbitrarily close to — Hence, there are models 
that start with a big bang and end with an asymptotic 
de Sitter expansion, but have an intermediate phase 
where d is very close to zero, so a(t) remains close to 
a for a very long time (solid line in Fig. |H]) and there 
is a phase in which the universe "coasts" Einstein's 
static model. This possibility was also discovered by 
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(4.8) 



whereas Eq. (|4.4p . together with the condition V{a) = 
— gives 



ao 



l + 3w 
{l + w)A 



1/2 



(4.9) 



Hence, the value of A determines uniquely both the 
scale factor and the energy density in this model. 
For w = one recovers Einstein's static universe<^ 



Time dependence of the scale factor in coasting 



FIG 

models. Solid curve: 
Dashed line: V{a) > 



V{a) 



/2 (Lemaitre universe) 



LemaitrC)^ see Ref. y for a detailed discussion. Of 
course, other models exhibiting such a coasting can 
be found, all corresponding to phase curves close to 
the separatrix (see, e.g., the dashed line in Fig. |8l). 

In the case fc = 0, it is also possible for the max- 
imum value of V to be very close to (i.e., to E). 
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However, the only possibility for a static solution in 
this case is a = (which happens when the ratio A/B 
attains the degenerate value A/B — 0). Therefore, a 
quasi-static phase in which a is almost constant, like in 
Lemaitre coasting models, is possible only right after 
the big bang, when the universe is extremely dense. 
At the end of such a phase, the universe turns even- 
tually into a de Sitter model. 

B. -1s;tos;-1/3 

The cases —1/3 < w < —2/3, w — —2/3, and 
—2/3 < w < —1, although mathematically distinct, 
exhibit many similarities. They all correspond to a 
fluid of topological defects (strings or domain walls).— 
We do not discuss all the different possibilities, but 
present only the case with — l<w<— 2/3 and A < 0, 
which exhibits interesting behaviours. It is an easy 
exercise to draw the phase curves and discover the 
behaviour of a{t) for other values of the parameters. 

The potential energy and the corresponding phase 
curves are plotted in Fig. [31 Figure [TU] shows the time 
behaviour of a. For k — —1, there is always a big bang 
followed by a big crunch. Note that, contrary to what 
happened when w > —1/3, now d is finite at these 
events. For A; = (separatrix) there is either a phys- 
ically uninteresting singular "universe" with a{t) = 
(dot at the origin of the phase plane), or a solution 
where the big bang and big crunch take place infinitely 
far in time. The behaviour in the case k = 1 de- 
pends on the minimum value of the potential energy, 
attained at a given by Eq. (|4.3|) . Such minimum value 
is, of course, still given by Eq. (|4.4|) . If V{a) > -c^/2 
(dashed curve for the potential energy in Fig. [S]) , no 
solution exists. If V{a) = — c^/2, there is a stable 
static model with scale factor given by Eq. (|4.9|t — 
(dotted curves for V{a) and a{t), and dot with a ^ 
in the phase portrait). Finally, if V{a) < —(? jl (solid 
curve for the potential energy) the universe oscillates 
around the static model. 

Let us now consider the border cases w = —\ and 
w — —1/3. For a fluid with w = —\ the stress-energy- 
momentum tensor— 

Tab = {p + p)UaUb+pgab (4-10) 

reduces to Tab = —pgab- This circumstance is inter- 
esting for two reasons. First of all, the four-velocity 
of the fluid does not appear in Tab, so the observable 
properties of the fluid do not allow one to identify a 
preferred reference frame. In particular, this is the 
case for the quantum vacuum,— ^—i^ which must, of 
course, be Lorentz-invariant. Second, the condition 
V^Tab — 0, expressing the local conservation of energy 
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FIG. 9: Potential energy and phase portrait for — 1 < u) < 
—2/3, A < 0, and three values of A. With the exception 
of the dot at a ^ 0, which denotes the stable static model 
when V is represented by the dotted curve, only the phase 
curves for the largest value of A have been drawn (sohd 
curve V{a)), because the same qualitative behaviour will 
occur in the other cases. 

and momentum, implies that p is a constant. '^^ Insert- 
ing the stress-energy-momentum tensor Tab — — Pgab 
into Einstein's equations, one then sees that such a 
fluid behaves like a cosmological constant with value 
A* = 8ttGpo/c^. The evolution of the universe is then 
entirely controlled by the effective cosmological con- 
stant Ac = A -t- A* as in Eq. (13.51) , and the potential 
energy becomes 

V{a)^-^a\ (4.11) 
Similarly, for w — —1/3, corresponding to a universe 
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FIG. 10: Time dependence of the scale factor correspond- 
ing to the phase curves of Fig. [O] 



dominated by a network of cosmic strings^ the po- 
tential energy is 



Via) 



Ac2 
~6~ 



(4.12) 



We do not discuss the behaviours corresponding to 
these simple forms of V{a), which are well-known from 
elementary mechanics. 



w < 



This is the "phantom energy" regime,— i^i^ intro- 
duced as an explanation for the observed acceleration 
of the uni verse. -iS>^ Now 1 -I- 3w < —2 so, defining 
6 = -3(l + w)/2 > such that (l + 3u;) = -2(1 + 6), 
the first contribution to potential energy is 



(4.13) 



Regardless of the value of A, this term dominates the 
potential energy for large values of a. In that regime, 
the total energy E is also negligible, and one can write 
the asymptotic differential equation 



2 



whose solution is 

a{t)^ (±SV2A{t^^~t] 



-i/s 



(4.14) 



(4.15) 



where t^^ are constants of integration. Hence, at some 
finite time in the future (positive signs in Eq. I4.15P or 
in the past (negative signs), the scale factor becomes 
infinite. The former event is usually referred to as the 
"big rip",^ 



It is interesting to estimate t^^.. If to is a time 
at which phantom energy dominates over all other 
forms of matter (including, possibly, a cosmological 
constant), equation (|4.15p with the positive signs is 
correct with excellent approximation at time Iq, so 



ao 



-1/5 



Hence, 



"br 




1 f8TTGpo\ 
d\3c^ J 



-1/2 



(4.16) 



(4.17) 



where we have substituted A from the first equation 
in (|3.4I) . If the cosmological constant and the spatial 
curvature (proportional to the numerical parameter k) 
are negligible at time to, Eqs. (|2.2|) and ()4.17|) imply 
that 



tn 



1 flo 



6 aito)' 



(4.18) 



Hence, if w is not too close to —1 (so S is not too 
close to zero), t^^. — to is of the same order of mag- 
nitude of the quantity ao/a(to)- This is the inverse 
of the Hubble parameter, which is comparable, in the 
standard model of cosmology, with the time elapsed 
since the big bangi^ Choosing to of the order of the 
present time (at which the cosmic acceleration starts 
to be relevant), it then turns out that t^^ — to is of 
the same order of magnitude of the present age of the 
universe (or smaller, if w has a large negative value). 
Therefore, the big rip happens in a not-too far future 
(in reference [13 an estimate of about 20 billion years 
for tj^j. — to is given assuming w = —3/2). On the other 
hand, for w very close to —1, 6 is very close to zero, 
and the big rip is delayed to a remote futurci^ 

It may seem pointless to consider combinations of 
A and phantom energy, which in cosmological models 
of current interest both serve to explain cosmic ac- 
celeration at the present epoch. However, there are 
good independent reasons to believe that A ^ 0^^— 
and if one admits the possibility of a component with 
w < —1 one can envisage models in which at earlier 
times A dominates over phantom energy, which later 
takes over and controls the time evolution in the re- 
mote future. For this reason, we present all such cases, 
again with the obvious caveat that none of them ap- 
plies to the evolution of the universe in its entirety. 



1. A ^ 

The potential energy and the phase portrait are 
shown in Fig. [TTJ and the time evolution of the scale 
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V 



a 




FIG. 11: Potential energy and phase portrait for w < —1, 
A ^ 0. 



factor is plotted in Fig. [121 For fc = — 1 the universe 
starts with a big bang and ends with a big rip, while 
for fe = (separatrix) the big bang happens infinitely 
far in the past (of course, mathematically, also the 
time- reversed processes are possible) . For k = there 
is also the unphysical possibility of a static singular 
"universe" with a{t) = 0, represented by the dot at 
the origin of the phase plane. For k = 1 the universe 
starts with an "anti-big rip" , contracts to a minimum 
size in a finite time and then re-expands symmetri- 
cally. 



a 




FIG. 12: Time dependence of the scale factor correspond- 
ing to the phase curves of Fig. 1111 Note the presence 
of vertical asymptotes (corresponding to the ticks on the 
t-axis) associated with a big rip and its time-reversed ver- 
sion. 



2. A < 

The potential energy and the phase portrait are 
shown in Fig. [T31 The time evolution of the scale 
factor is plotted in Fig. [141 For k = 1 and k — 0, 
the universe starts with an "anti-big rip" , contracts 
to a minimum size, and then re-expands symmetri- 
cally until a big rip occurs. (Again, for k ^ there is 
also a static singular solution a(t) = 0, corresponding 
to the dot at the origin of the phase plane.) For the 
value (|4.3p of the scale factor, the potential energy 
has a maximum, given by Eq. (14.41) .'^ Hence, there 
are three possibilities for k ~ —1, which lead to qual- 
itatively different behaviours. 

If V{a) > c^/2, there are two disconnected ranges 
for the possible values of a, corresponding either to a 
big bang followed by a big crunch, or to a model in 
which a contracts from infinity to a minimum size in 
a finite time and then re-expands symmetrically, as in 
the cases with k — 1 and fc = 0. 

If V{a) = we have a separatrix (dotted curves 
in Figs. [131 and [T4l) . There is an unstable static so- 
lution a{t) = a, analogous to Einstein's static model 
of Sec. IIVA3I represented by the dot with a ^ in 
the phase portrait. There are also a solution that 
starts with a big bang and approaches the static model 
asymptotically for t — >■ -l-oo, and one that coasts the 
static model for an infinite time and then ends with a 
big rip (as well as their time- reversed versions). 

Finally, if V{a) < c^/2 (dashed curves), the uni- 
verse starts with a big bang and ends with a big rip 
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a 




a 




FIG. 13: Potential energy and phase portrait for w < —1, 
A < 0, and three values of A. The phase curves have been 
drawn for all the three values of k only for the smaller value 
of A (solid curve V{a), solid phase curves). For the other 
potential energies, only the phase curves corresponding to 
fc = — 1 are plotted. 



(or the time- reversed process). Interestingly, during 
this evolution a changes sign from negative to posi- 
tive. If V(a) is very close to the critical value 
there is a phase in which a{t) is approximately equal 
to the value a and increases very slowly, as in the 
Lemaitre coasting models of Fig. HI 




FIG. 14: Time dependence of the scale factor correspond- 
ing to the phase curves of Fig. 1131 Note the presence of 
vertical asymptotes where a — >■ -|-oo at a finite time. 



V. FINAL REMARKS 

Using standard techniques of qualitative study, fa- 
miliar from courses in classical mechanics, we have in- 
vestigated the time evolution of Friedmann-Lemaitre- 
Robertson- Walker models with a perfect fluid and a 
cosmological constant. The analysis is elementary and 
does not require one to solve differential equations (ex- 
cept in some very simple cases, if one also wants to 
know the detailed time behaviour in some asymptotic 
regimes). 

There are several ways in which our discussion could 
be expanded. Although we have considered only cases 
in which the parameter w in the equation of state (11.11) 
is flxed for all times, one can easily construct "patch- 
work" models where the cosmic evolution is domi- 
nated by different kinds of matter at different eras. 
Also, one can combine the behaviour found for a with 
Eq. (|2.6I) in order to find the time evolution of the 
energy density p. Furthermore, the qualitative anal- 
ysis can be straightforwardly extended to a generic 
barotropic equation of state (|2.4[) . Indeed, integrating 
Eq. (ESI) gives 



dp 



P + P{p) 



(5.1) 



from which one finds p as a function of a. This can 
then be replaced into the Friedmann equation (|2.2|) to 
obtain the potential energy for the equivalent mechan- 
ical system. For example, a generalised Chaplygin gas 
with equation of state 



(5.2) 
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with C > and < a ^ 1, has sometimes been con- 
sidered as an alternative to phantom energy, in order 
to explain the cosmic acceleration .^^'^^ This form of 
p{p) leads to 

with D = (Po^^ — C) aQ^"^^\ and thus to a contribu- 
tion 

«^ (5-4) 

to the potential energy V{a) in Eq. p.ip . It is again a 
simple exercise to find the asymptotic behaviours and 
to perform the qualitative analysis. 

As an alternative to qualitative analysis, one could 
integrate numerically the Friedmann equation (|2.7p 
to produce the plot for a(t). Although one might well 
prefer such a strategy in professional research work, 
we believe our approach to be pedagogically more 
valuable. With the help of a phase portrait, a stu- 
dent knows at a glance the behaviour of a cosmological 
model under all possible initial conditions, recovering 
also cases that are likely to escape when one performs 
only a few numerical experiments. (This happens, 
e.g., close to a separatrix; for example, it would be 
hard to hit numerically on the Lemaitre models of 
Fig. |S1 unless one knows already what to look for.) 
In order to produce good-quality figures, all the plots 
in this article have indeed been obtained numerically. 
However, it is important to realise that they could 



be reproduced with all their essential features, using 
only pencil and paper and a few techniques known 
from elementary calculus, so the technique presented 
here allows students to quickly and easily grasp the 
main features of a cosmological model. 

Of course, not all the combinations of w and A 
are equally important for an understanding of the 
standard model of cosmology. Only the content of 
Sec. IIV Al (w > —1/3) is mandatory for an introduc- 
tory course. However, as we have seen, the generality 
of the method allows an instructor to easily present 
some conclusions pertaining to other, more exotic, val- 
ues of w whenever that might be useful. 

The material presented in this paper can be used 
not only in a course on cosmology, but also as an un- 
usual illustration of the qualitative techniques during 
lectures on classical dynamics, thus allowing an inter- 
esting interplay between concepts from different areas 
of physics. Although the analysis follows the tradi- 
tional scheme, the cosmological interpretation of the 
variable a, and the corresponding description in terms 
of the behaviour of models of the universe, should pro- 
vide additional motivations for the students, who of- 
ten find these subjects more fascinating and attractive 
than the point particles of classical mechanics. 
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